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“The aeronautical and space activities of the United States shall be 
conducted so as to contribute... to the expamsion of human knowl- 
edge of phenomena in the aimosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 
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TECHNICAL REPORTS: Scientific and technical information considered et 
important, complete, and a lasting contribution to existing knowledge. H ' 
TECHNICAL NOTES: Information less bread in scope bur nevertheless of : ‘4 i 
importance as a contribution to existing knowledge. -&. 


TECHNICAL MEMORANDUMS: Information receiving limited distribu- ote 
tion because of preliminary data, security classification, or other reasons. wf 


CONTRACTOR REPORTS: Scientific and technical isformation generated 
under a NASA contract or grant and considered an impertant contribution to 
existing knowledge. 


TECHNICAL TRANSLATIONS: Information published in a foreiga {- 
language considered to merit NASA distribution in English. b 


SPECIAL PUBLICATIONS: Information derived from cr of value to NASA 
activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 


FECHNOLCGY UTILIZATION PUBLICATIONS: Information on tech: 


Be 
x 
wa “2, 
Poe, Haden retecbentsgiO ye 


gee 
oats. 


Bain 


nology used by NASA that may be of particular interest in commercial and other } 

non-atrospace applications. Publications include Tech Briefs, Technology 3 i 
Usilization Reports and Notes, and Technology Surveys. < 
VE: 
a@ i 
Details on the availability of these publications may be objained from: y 4 : 
ee 
SCIENTIFIC AND TECHNICAL INFORMATION DIVISION H 
4 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION qi ; 

Washington, D.C, 20546 i 
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Preface 


- A program sponsored jointly by the National Aeronautics and Space 
Administration and the United States Air Force has been in progress 
under the direction of Dr. Vette for the purpose of defining a model 
radiation environment of the Earth. In Volumes I and II of NASA 
SP-3024, the environment was given for the lower altitude region where 
trapping is relatively stable and changes in radiation flux occur only 
slowly and are generally small, In Volume MI, the electron environ- 
ment was described which is found at 19,300 n. miles, the altitude of 
the Earth synchronous orbit. The present report contains a mode} of 
the low energy proton flux for energies below 4 MeV. 


Since low energy protons are stopped very easily, their contribution 
to radiation damage can be neglected in most cases. Over a iarge re- 
gion of space, however, their flux is intense, thus they may be the prime 
cause of deterioration of sensitive surfaces. These effects may be 
particularly important to temperature control surfaces, optical windows, 
or the first elementin a counterielescope for studying nuclear radiation. 


This compilation would not have been possible without the assist- 
ance given by the investigator who performed the original measure- 
ments. All users of this model environment will greatly appreciate 
these efforts. Clearly, no model can be better than the data available 
for deriving it. 


A. W. Schardt 


Physics and Astronomy Programs 
NASA Headquarters 
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INTRODUCTION 


This report is the continuation of a series of model environmerts of the particle radiations 
trapped in the geomagnetic field. Preceding environments have described protons of energies 
greater than 4 MeV trapped in the inner zone and electrons trapped in the inner and outer zones 
and at synchronous altitudes (References 1, 2, 3). The four previous proton environments are dis- 


tinguished hy .ne energy range over which each is applicable. The omnidirectional, integral fluxes 
are presented in the form 


JOE: BL) = J(>E,; BL) eo (EOE1)/E0B.t) 


where E, is the reference energy fer each environment. The spectral parameter E, (B, L) is dif- 
ferent in the various environments, or equivalently, EF, (B, L) is a discontinuous function of energy. 
The criterion for determining the energy range over which an environment may be applied haa 
been the existence of an exponertial parameter E, which represents the data over that range with 
the desired accuracy. It was found necessary to construct previous environments over the energy 
ranges 4 - 15 MeV, 15 - 30 MeV, 30 - 50 MeV, and >50 MeV. 


The present environment, AP5, describe3 protons of energies below 4.0 MeV trapped between 
L = 1.2 and L = 6.6. The general approach is the same as in the preceding proton environments, 
but there are two main differences. Again the energy range of validity is that over which a simple 
spectral representation suffices, but we inquire into whether a power law cependence or an ex~ 
ponential dependence is more representative of the data below 4 MeV. (Neither representation will 
be more accurate ihan the approximate factor of two uncertainty associated with the temporal 
fluctuations in the data.) It is found that, over the complete region of 8-£ space of interest, neither 
representation provides a significantly better fit to the data than the other. Consequently the ex- 
ponential representation is chosen for the sake of continuity with the preceding environments. 


The second difference between this and previcus environments is that in this case there is not 
a unique energy which may be assigned as the lower energy limit. As will be discussed more fully, 
there is an L-dependent peak in the differential energy spectrum which leads to an L-dependent 
lower energy limit of the model. This lower limit may be taken as 0.1 MeV for L > 3.5, and grad- 
ually increases as L decreases. (See Figure 3.) 


The distribution function J(>0.4 MeV; B, L) and spectral function E, (B, L) are presented in 
tahular form. Many graphical comparisons of the model with the data are given. Orbital integra- 


tion tables are presented for circular orbits up to 18,000 n. miles and for 15 energies and 14 energy 
bands. 
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TEMPORAL VARIATIONS 


In Table 1 is presented alist of the data employed in the construction of this environment, All 
the experiments involvedthe measurement of unidirectional proton fluxes, Scintillation deteciars 
were utilized on the P-11 and the three Exnlorer experiments, while solid state detectors were 
used en the Relay 1 and Injun 1 experiments. During the period over which the data coverage ex- 
tends (July 1961 - April 1965), temporal variations of up to a factor of two occur in the measured 
fluxes, at least for L< 4.5. The growth and decay factors involved in these changes are position- 
and energy-dependent. That the temporal changes are fluctuations and not monotonic processes 
as in the case of Starfish electron decay is indicated by the data. Between December 1962 and 
May 1963, Fillius and McIlwain observe a slight decay in the fluxes of protons with energies 
greater than 1.1 MeV (Reference 5). Two years later White, and Davis and Williamson observed 
an increase in the fluxes of these particles and a decrease in the fluxes of lower energy protons 
(References 9, 11). Detailed discussions of the observed temporal dependence of low energy pro- 
ton fluxes are presented in References 6, 9 and 11, 


As larger L values are considered, fluctuation arnplitudes also become larger. Davis and 
Williamson have reported on rapid (ten minutes) order of magnitude changes of proton fluxes at 
L = 4.9 associated with the 17 April 1965 magnetic storm (Reference 9). Also, the 1962 Davis- 
Williamson data (flux greater than 0.5 MeV) are in apparent disagreement with the corresponding 
data obtained two years later on Mariner IV (Reference 12). The lattcr is about two orders of mag- 
nitude smaller than the former fo~ L 26.0. 


This model environment is presented with no explicit time dependence, despite the foregoing 
discussion. At L <4.5 there are not yet sufficient data to determine the long term (say, eleven 
years) temporal behavior of the proton fluxes, but there are sufficient data to establish the fluc- 
tuating character of the temporal variations and to establish the factor of two as an approximate 
upper limit to the fluctuation amplitude. Thus the fluxes predicted by the model should be within 
a factor of two cf any of the observed fluxes. Graphs to be cons.dered later will demonstrate the 
validity of this statement. 


At L 2 4.5, we have chosen to utilize the 1962 Davis-Williamson data near the equator rather 
than the Krimigis~Armstrong Mariner IV data. The former are more extensive and indicate higher 
flux values, which are perhaps more appropriate for applications purposes. At large B values the 
Mihalov-White data are utilized, and interpolated flux values are assigned to intermediate B values. 
When sufficient low energy proton data are available in this high L region, a time-statistical anal- 
ysis of the proton fluxes (comparable to that performed on synchronous electrons in AE3) may be 
warranted. It is in the region L 25.0 that the ability of the B-L coordinate system to order the data 
begins to fail due to the neglect of magnetospheric distortion. These factors should be borne in 
mind in using the environment at very large L values. 
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jAPS ENERGY LIMITS 


Peaked differential energy spectra of low energy proton fluxes have been observed by Mihalov 

and White (Reference 8) and by Davis and Williamson (Reference 9). As L increases, the energy 

at which the peak occurs decreases. This behavior is also theoretically predicted by diffusion 

theory in which the first two adiabatic invariants, but not the third, are conserved in particle mo- 

tion (Reference 13), Shown in Figure 1 is a plot of the energies at which this peak was observed at 

Se different L values. At a given L value, the peak occurs at larger energy as B decreases (i. e., as 
the equator is approached). 


It is desired to include in this environment only that energy range in which the data can be 
‘represented by a simple spectral function (power law or exponential) to within the desired degree 
of accuracy. Because the accuracy of the model has been limited to a factor of two by the time 
fluctuations in the data, the lower energy limit of this model may be extended to energies some- 
what (but not arbitrarily) lower than the peak energy of the differential spe~‘rum. 


| 
i 
} 
In Table 1 are presented the L values over which the various sets of data have been used. 
Curves to be presented later show that, at any given L value, the data do agree with the model to 
within a factor of two. The lower energy limit to be assigned to the model at’ given L is lower 
than that of the data utilized at that L value, as listed in Table 1. To illustrate this, we present in 
Figure 2 the environment flux at the equator, where the inclusion of very low energy data is least 
likely to match the environment. A comparison of the 1962 and 1965 Davis-Williainson fluxes of 
protons (E > 0.1 MeV) with the environment is presented. It is clear that, except near L = 4.5, 
the E>'0.1 MeV data agree with the model beyond L * 3.5. We take 0.1 MeV to be the low energy 
limit of the environment beyond L = 3.5. Data on protons of energies lower than 0.1 MeV are not 
yet available. 
| 
i 


In the region L< 2.0 only the Relay 1 data and the Injun I data were available at the time the 
environment was constructed. Neither of these involves proton energies less than 1.0 MeV. The 
data of Fillius and McIlwain suggest that the differential spectral peak occurs at energies above 
1.1 MeV in this region (Reference 5), as may have been expected by an extension tc L < 2.0 of the 

curve determined by the Davis equatorial points of Figure 1. The nested character of the energy 
intervals in the Fiullius~Mcilwain experiment renders impossible a more accurate determination 
of peak energy. Despite this, a single spectral representation is assumed as appropriate for this 
data, and ag shown in Figure 27 that even with this approximation. the model is representative of 
the data to within a factor of two. Below L = 1.5 there are no comprehensive data available in the 
energy range of interest; as such the flux values attributed to the model for L <1.5 are extrapolated 
from the data available at L >1.5. 


Presented in Figure 3 is a curve of the L-dependent energies which may be assigned as the 
lower énergy limit of the model while yielding an accuracy of a factor of two for the 
environment. 


*2 JE; B,L), The function J (>E; 8, L) is presented as the product of a distribution function of fluxes 


PROTON ENERGY SPECTRUM 
The environment is presented in terms of integral, omnidirectional ‘luxes J (>E; B, L, t) 


above the reference energy, E,, and a spectral function, N(>E; B, L): 
JOE; B,L) = J(>Ey; B,L)NOE;B.L) . 
It is assumed that the spectral function is of a simple form. Both a power law representation, 
N(PE; B, L) = (B/E,)?@4) , 


and an exponential representation, 


NODE; B,L) = @ PED PORE 


were examined. That representation which better fit al’ the da --er the entire B-L range of the 
environment was to be chosen. However, as wili be shown, the spz-.d of the data does not clearly 
dictate the choice cf one representation over the other, even though some of the data favors an 
exponential spectrum (Reference 4) while other data favors a power law spectrum (Reference 8). 
The temporal fluctuations in the data indicate that the consideration of a spectral function of any 
greater complexity is unwarranted. 


The manner of determining the better representation and best parameter (p, E,) ata given 
point in B-L space is as follows. Let J; (E,2~-E,,; 8, L) be the it* observed flux value at the B, L 
point, representing a flux of particles lying in the energy interval E,z to E,,, (where Ej, may be 
infinity). If an exponential spectrum is assumed, such that 


JOE: B,L) = J(PE,3B, L) oe EEE o(BL) 


then the flux at the reference energy, E., corresponding to the i‘* data value and j'" possible E, 
is given by 


~E,/E,, (BL 
ooE1/Fo3(B-L) 


i ~E£/Fo; .B.L) ng Fiv/®o} (2, » 


J;; PEB.L) = J; (E.2-8,,3 BL) 
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Alternatively, if a power law spectrum is assumed, such that J(?E; B, L) = J(> Ey B,L)- (E/E,) 78 , 
then the flux at the reference energy corresponding to the it" data value and j'" possible pis 
) es ha 
>E,3B, = EQ°E 8 be 
Ji; ( E33, L J; ( i27E,,:B L) errr -E. as) 


If we assume that there are n observed data values at the B, L point of interest, then the flux 


chosen at the reference energy (i.e., the distribution function) is the geometric mean of the J 
that is 


ij? 


si/n 
J, PEiBL) = rr Ji; (7E.sB, | . 
i=l 


This, of course, is equivalent to 


1 n 
log]; (E;:BeL) = 5 ye logJ,; (E,;B,L) . 
i=1 


We determine the better representation and best spectral function as those which minimize 
the logarithmic root mean square deviation of the fluxes at the reference energy, 


rs 1/2 
a oe 2amt (toe 5, (PE,; B. L) - log J,, (°E,: B, L))" . 
i=1 : 


Minimizing o, for each representation yields the best spectral parameter for each spectral repre- 
sentation. The better representation is that one giving the smaller minimum ¢. Thus the spectral 
representation and parameter which best describes all the data are selected at each of a sufficiently 


fixe grid of points in B-L space. A computer program was generated by A. B. Lucero to perform 
the preceding analysis. 


After the best E, and best p are chosen at each B-L point, curves of E, versus B and of p versus 
B are drawn at fixed L values. Because the various sets of data overlap only partially in B-L space, 
a certain amount of smoothing of these curves is necessary. Presented in Figures 4 ~ 21 are the 
E, and p values at selected L values. The dots are obtained from the pointwise consideration of the 
data, while the solid lines are the smoothed values. 


It is possible to determine spectral parameters by a consideration of data obtained un one ex- 
periment only. Such parameters have been determined and are represented by the various symbols 
shown in rigures 4 - 21. The differences among these values may be jue in part to the differing 
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energies or observation times involved. They are included to give the reader an estimate of the 
disagreement among the various experiments. 


The procedure actually followed was to consider both the power law and exponential spectral 
representations, with the smoothed values of E, and p used to determine distribution fuactions. 
The logarithmic root mean square deviations (c) were again obtained for each representation, and 
the values of o were averaged over B for fixed L (i.e., along a line of force). The results of this 
calculation are shown in Figure 22. A calculation of 


0 
I= os (av)i , 
ist 


where (c,,), is the average value of cat L = 2+ .2i, yields the result that I = .78 for the ex- 
ponential representation, while 1 = .80 for the power law representation. It is clear that over 

the complete range of interest, neither representation gives results significantly better than the 
other. As such, we have elected to present the details of this environment using an exponential 
spectrum, because preceding environments have utilized exponential spectra. The values of E, which 
constitute the spectral part of this environment are presented in Figure 23 and in Table I. It is 
clear that the spectrum softens as L increases and that, at a fixed L value, the spectrum generally 
softens as B increases. 


DISTRIBUTION FUNCTION 


After the distribution function J(>E, = 0.4 MeV; B, L) is obtained as described previously, a 
log B-L plot of isoflux contours is drawr and smoothed. The result is a distribution function 
smoothed in both Band L. This log B-L plot is shown in Figure 24, while the distribution function 
is given in tabular form in Table 2. Figure 25 shows the distribution function in the more physically 
intuitive R-\ representation. It is apparent that there is a sharp decrease in the flux levels as the 
atmosphere is approached along a given magnetic field line. In Figure 26 the equatorial fluxes of 
protons with energies above 0.4 MeV and 1.0 MeV are presented. Maxima in these curves occur 
near L = 3.0 and L = 2.8 respectively. The 2> 0.4 MeV fluxes are based in part on data taken 
near 0.4 MeV for L > 2.0 but are based on extrapolations from data taken on E > 1.0 MeV proton 
fluxes for L < 2.0. This discontinuity in the data leads to the slight irregularity near L = 2.0 in 
the E >0.4 MeV curv? of Figure 26. Data on protons with energies above 1.0 MeV were incorporated 
into the model for L 21.5, and this is reflected in the smoothness of the E >1.0 MeV curve of Fig- 
ure 26. 


Figures 27 - 37 show the distribution function at selected L values, along with the values of 
the experimentally observed fluxes as reduced to the reference energy by means of the model ex- 
ponential spectrum. It may be observed that the environment agrees with the data to within a fac- 
tor of two or better, nearly everywhere. Shown in Figures 38 ~- 42 are the distribution function and 
corresponding experimental fluxes obtained at integer L values using a power law spectral 


2 NRE AEA TRESS RINE ACT EAE OEE SERRE RRNA I, APRA Fe 


J 


: 
| 
i 
| 
: 
i 
4 
4 
| 
5 
ir 
i 
1 
' 
: 
i 
Hu 


wry Rae Ke Ufa Gaba BNC Rinse soe comin « Ce ee me PSM Wes embe NAF AH Ue ERR oak has TE BSEA oH GLAM a HEP RE RAE re ene A eee 
. ae 


ney rs 
[saat 


| 
| 
| 


. { see ronan ine seh an ea mt nme ee ne et ae te we 


sete cenetiee 


2-5 ~~ a. 


a 


. SS eer er eee 
Fi 


representation. Comparison of Figures 27 ~ 3% with Figures 38 - 42 again illustrates the fact that 
neither spectral representation is significantly preferred by the data. 


Fiuxes of (= > 0.4 MeV) protons circulating perpendicular to the magnetic field direction have 
been obtained by conversion from the APS omnidirectional fluxes. These unidirectional fluxes are 
presented in graphical form in Figure 43 and in tabular form in Table 2. 


ORBITAL INTEGRATIONS 


Fluxes accumulated in various satellite orbits have been obtained by orbital integrations dis- 
cussed in detail in Reference 1. Orbital integrations have been performed for circular orbits up 
to 18,000 n. miles in altitude and at inclinations of 0°, 30°, 60° and 90°. Fluxes accumulated in 14 
energy bands-and the integral fluxes above the 15 associated energy levels have been calculated on 
a computer. These results are presented in Tab'e 3, with the flux values in units of protons/cm?- 
day. The time interval between successive orbit points and the total running time are shown; 
however, the total accumulated fluxes have been divided by the running time, in days, to obtain the 
average flux per day. (The symbol * is to be read "greater than.") Figure 44 illustrates the ac- 
cumulated daity fluxes of protons with energies above 0.4 MeV. 
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DISCUSSION 


we 


A model environment of the fluxes of protons having energies below 4.0 MeV has been pre- 
sented. Because of the peaked character of the differential energy spectrum, it was necessary to 
assign an L-dependent lower energy limit to the model. This is shown in Figure 3. The data 
coverage was best in the region 2.4 <L<4.0, and only limited in the regions 1.2 <L<2.0 and 
5.2<L<6.8. Because of the temporal fluctuations in the data, the accuracy of the model is given 
approximately by a factor of two over much of the region of B, L space. Larger amplitude temporal 
variations and a relative shortage of data render the environment somewhat less accurate at L > 5.0. 
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It was determined that over the energy and B-L ranges of interest, neither the exponential nor 
the power law spectral dependence was significantly more representative of the data than the other. 
: An exponential representation was chosen because previous proton environments have utilized ex- 
rd 


ee 


ponential spectra. Individual experiments (e.g. Reference 8) have indicated that a ;ower law rep- 
resentation may permit the construction of individual environments valid over greater energy 


ranges than has been possible using exponential spectra. This will be investigated after all the 
proton environments have been updated. 
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Figure 1—Locations of peak in differential energy 
epectrum as observed by the indicated experimenters. 
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Figure 3~Recommended low energy limit of 
proton map APS. 
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Figure 5—Comperisen of spectral parameter Ey 
with data at L = 2.5. 


10 


10 14 (8 2.2 26 30 34 38 42 46 50 


<) 


% 


3, 


OMNIDIRECTIONAL FLUX (E204 Mev) 
AT THE GEOMAGNETIC EQUATOR 
APs 
© DAVIS, EXPLORERS 12, 4, £>00 8¥ 
X_ DAVIS.EXPLORER 26,01Me¥ 


3 


CRINDIRECTIONAL FLUX > 04 MeV (protans Aan -29¢) 


9 8 26 34 42 30 38 66 
LIEARTH RADII) 


Figure 2—Comporison of AP5 equatorial flux with fluxes 
of E > 0.1 MeV protons observed in 1962 and 1965; the 
AP5 spectru:n has been used. 
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tained from one experiment. Solid line is E(B) used 
in AP5. 
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Figure 6—Compaorison of spectral poremeter Eq 
with dota at L = 3.0. 
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Figure 7~—Comparison of spectral parameter Eo 
with date at L = 3.5. 
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Figure 11—Comporison of spectral parameter Ep 
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Figure 8—Comparison of spectral parameter E, 
with data at L = 4.0. 
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Figure 10—Comparison of spectral parameter Ey 


with data at L = 5.0. 
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Figure 12—Comporison of spectral parameter Ep 
with data ot L = 6.0. 
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Figure 13—Comporison of power law parameter p 
with data at L = 2.0. 


PORE PEAR ERG 
L930 


@ ALL DATA 
O FILLES Mc RUN, RELAY | 


O WHALE /-WHITE, PIIE<LS Me) 
° MIKALOY- WHITE, PHIE>LS Mev) 


Figure 15—Comparison of power law poramete: p 
with data at L = 3.0. 
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Figure 17—Comparison of power law poremeter p 
with data at L = 4.0. 
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Figure 14—Comporison of power law porometer p 
with date at L = 2.5. 
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Figure 146—Comparison of power law porameter ¢ 
with data at L = 3.5. 
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Figure 18—Comparison of power law parameter p 
with data at L = 4.5. 
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Figure 19—Comparison of power law parameter p 
with dota at L = 5.0. 


POWER L&Y PARAMETER p 


& DAVIS - WELIAMSON, EXPLOCERS 12,14 
O MIHALOV- SHITE, PU 


Figure 21—Comparison of power law parameter p 
with data at L = 6.0. 
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Figure 20—Comparison of power law porameter p 
with data ot L = 5.5. 
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for fixed L, using the two spectral representations. See 
text for definition of o. 
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Figure 23—Spectral parameter E, used in the proton environment AP5. 
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Figure 24—Log B~L flux mop of the APS environment. 
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Figure 26—Omnidirectional proton flux at 
the geomagnetic equator. 
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Figure 25—R-A flux map of the AP5 environment. 
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Figure 27—Comparison of proton map APS 
with data at L = 1.6. 
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with date at L = 2.0. with data at L = 3.5. 
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Figure 32—Comparison of proton map AP5 
with data at L = 4.0. 
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Figure 34—Comparison of proton map AP5 
with dota at L = 5.0. 


OMNIDIRECTIONEL FLUX > 0.4 MeV (prutons /om? -sec) 


B(gauss) 


Figure %3—Comparison of proton map AP5 
with data ot L = 4.5. 
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Figure 35--Comparison of proton map AP5 
with data ot L = 5.5. 
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GESTRIBUTION FUNCTION FOR 
POWER LAK SPECTRUM 
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Figure 40—Comporison of data with distri- 
bution function obtained at L = 4.0 using 
power law energy spectrum. 
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Figure 42—Comparisor of data with distri- 
bution function obtained at L = 6.0 usin.a 
power low energy spectrum. 
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Figure 41--Comparison of data with distri- 
bution function obtained at L = 5.0 using 
power law energy spectrum. 
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Figure 43—-Unidirectional flux perpendiculer 
to magnetic field obtained fron: AP5 distri- 
bution function. 
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